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Multiple  Impulse  Solutions  to  McKean’s  Caricature 
of  the  Nerve  Equation.  II.  Stability 
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WEI-PING  WANG 

The  University  of  North  Carolina  at  Chapel  Hill 
Abstract 


We  studyMcKean’s  caricature  of  a  nerve  conduction  equation 


(1) 


dip  d2u  .  .  . 

it~  j?'u  +  H(u~°)  +  v’ 

-  j,,  .  . 


0  <  a  <  1, 


dv  . 

jj=bu-cv. 


b  <  0.  c  >  0, 


where  H  is  the  Heaviside  function.  It  is  proved  that  an  n-ple  impulse  solution  resembling  the 
superposition  of  n  unstable  solitary  impulses  has  at  most  2n  -  1,  and  at  least  n.  unstable  modes: 
exactly  n  unstable  modes  corresponding  to  the  amplitudes  and  the  rest  of  them  corresponding  to  the 
spacings.  The  n  amplitude  modes  always  exist.  W«.-prove  also  that  for  an  n-ple  impulse  solution 
resembling  the  superposition  of  n  stable  solitary  impulsesf  there  are  at  most  n  -  1  unstable  modes 
and  all  of  them  are  of  spacing  type.  /  .  (  I  ’ 

/  -  ,r  '-'J 

■  / 

Introduction 

In  the  first  part  of  this  series  [15],  we  proved  that  for  fixed  a,  h,  c  (in  a 
reasonable  range)  and  arbitrary  positive  integer  n,  the  system  (1)  has  countably 
many  multiple  impulse  solutions  consisting  of  n  widely  spaced  solitary  pulses. 

The  issue  of  the  stability  of  multiple  impulse  solutions  is  particularly  interest¬ 
ing.  At  some  initial  time,  an  impulse  is  subjected  to  a  perturbation  which  is 
bounded  spatially.  If  the  system  responds  by  causing  the  perturbation  to  grow  in 
time,  thereby  changing  the  form  of  the  impulse,  the  impulse  is  unstable.  If  the 
perturbation  decays  leaving  the  impulse  traveling  with  the  same  form,  it  is  stable. 

Stability  of  solitary  impulses  of  McKean’s  caricature  is  clear.  Rinzel  and 
Keller  [14]  demonstrated  by  computing  a  growth  rate  of  the  instability  that  the 
slow  solitary  impulse  is  unstable,  and  conjectured  that  the  fast  solitary  impulse  is 
stable.  Feroe  [3]  completed  Rinzel  and  Keller’s  work  by  computing  a  winding 
number.  [3]  proved  numerically  that  the  fast  solitary  impulse  is  stable  and  the 
slow  solitary  impulse  has  exactly  one  unstable  mode. 

Evans  [1]  studied  the  stability  of  multiple  impulse  solutions  for  the  general 
Hodgkin-Huxley  system.  He  found  a  complex  analytic  function  D(X)  defined  in 
a  domain  containing  the  right  half  of  the  complex  \-plane.  It  is  proved  in  [1]  that 
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(a)  Slow  solitary  impulse 
Winding  number  -  1. 


(b)  Fast  solitary  impulse 
Winding  number  -  0 


(c)  Slow  double  impulse. 
Winding  number  «  3 


(d)  Fast  double  impulse 
Winding  number  *  0 


(e) 


Slow  double  impulse 
Winding  number  =  2. 


(f)  Fast  double  impulse 
Winding  number  =  1 


Figure  1.  The  image  of  the  imaginary  axis  under  the  mapping  £>. 


the  winding  number  of  the  image  of  the  imaginary  axis  under  the  mapping  D  is 
the  same  as  the  number  of  unstable  modes.  Since  the  function  D(\)  is  extremely 
complicated,  the  winding  number  can  only  be  found  by  computer.  Feroe  applied 
Evan’s  winding  number  to  McKean’s  caricature  with  solitary  impulses  in  [3]  and 
a  few  double  impulses  in  [5],  The  image  of  the  function  D  for  n-ple  impulses, 
with  n  ^  3,  would  be  too  complicated  to  count.  However,  it  is  interesting  to  look 
at  Feroe  [3]  and  [5]’s  winding  number  picture,  with  n  <;  2,  for  McKean’s 
caricature  (see  Figure  l)1: 


Most  of  the  figures  were  kindly  supplied  by  Dr.  John  A.  Feroe,  with  the  permission  of  SIAM. 
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Feroe  [5]  also  sketched  the  pictures  of  the  unstable  modes  for  double  impulses 
and  conjectured  that  there  are  only  two  kinds  of  instabilities:  spacing  type  and 
amplitude  type. 

In  this  paper  we  give  an  analytical  answer  to  both  the  number  and  the  type  of 
the  unstable  modes,  of  multiple  impulse  solutions  consisting  of  any  finite  number 
of  pulses,  of  McKean’s  caricature.  The  results  are  summarized  in  the  following 
diagrams: 


I  The  number  of  unstable  modes  of  a 
\  “superposition”  of  n  slow  solitary  waves 

/  \ 


exactly  n  modes  of 

at  most  n  -  1  modes  of 

amplitude  type 

spacing  type 

j  The  number  of  unstable  modes  of  a 
\  “superposition”  of  n  fast  solitary  waves 

J 

they  are  all  of  spacing  type 


The  author  believes  that  the  upper  and  lower  bounds  provided  above  are 
crucial.  Proving  this  fact  is  an  interesting  problem  for  future  research.  As  a 
special  case,  which  is  of  most  importance,  when  the  lower  bound  zero  of  the 
second  inequality  is  reached,  one  obtains  a  stable  n-ple  impulse.  An  even  more 
important  open  problem  is  to  provide  criteria  that  determine  when  an  n-ple 
impulse  traveling  wave  is  actually  stable. 

All  of  the  results  mentioned  above  are  about  infinitesimal  stability.  It  would 
be  very  interesting  to  look  at  the  real  evolution  in  function  space  of  the  full 
nonlinear  equation.  McKean-Moll  [13]  studied  this  problem  for  a  “small” 
function  space  and  a  special  range  of  a,  b ,  c. 


1.  The  Variational  Equation 

We  introduce  the  traveling  coordinate  frame  (z,  t),  with  z  =  x  +  kt,  in  which 
system  (1)  takes  the  form 


du  _  d2u  .  du 
dt  dz2  k  dz 


-  u  +  H(u  -  a)  +  v. 


dv  _ 
It  ~ 


+  bu  -  cv. 


(2) 
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For  any  positive  integer  n,  the  n-ple  impulse  (uk,vk)  is  a  /-independent  solution 
of  this  system.  To  study  its  stability  we  consider  the  variational  equation 

(3) 

Since  uk(0)  =  uk(z2)  =  •••  =  uk(z2n_1)  =  a,  the  compound  6-function  is  ex¬ 
plicitly 

(4)  Huk-a)  =  2n£ 

We  now  look  for  solutions  to  the  system  (3)  having  the  form 

(5)  U(z,  t)  =  ex‘U(z),  V(z,t)  =  ex'V(z). 

It  follows  from  (3)  and  (4)  that  (V,  U,  V)  must  satisfy  the  ordinary  differential 
equation 


'U' 

r* 

1  +  X 

(6) 

u 

=  i 

0 

V 

L° 

b/k 

for  z  *  z,, 

i  =  0,1,2,- 

••,2 n  - 

1,  and  1 

(7) 

1/1 

£=- 

t/U) 
K  (*,-)! 

f  =  0,1,---,  2n  —  1 . 


The  H-ple  impulse  uk  is  unstable  if  the  system  (6)  and  (7)  has  a  bounded 
solution  with  dl*X  >  0.  This  solution  is  an  unstable  mode  with  growth  parameter 
X.  On  the  other  hand,  uk  is  stable  if  the  system  (6)  and  (7)  has  no  bounded 
solutions  with  die  X  >  0.  For  \  =  0,  there  is  always  a  solution  to  the  system  (6) 
and  (7):  (U\  U ,  V)  =  (uk,  u'k,  v'k).  This  is  because  any  translate  of  uk  is  also  a 
solution  to  the  system  (2). 


2.  The  Main  Results 

Our  rest  is  based  on  Rinzel-Keller  [14]’s  numerical  evidence,  i.e.,  we  assume 
the  a-k  curve  and  a-X  curve  there  (see  Figures  2  and  3).  Their  computation  was 
carried  out  for  the  case  c  =  0,  but  the  results  remain  true  when  c  is  small  by 
continuity.  The  a-\  curve  corresponds  to  the  lower  branch  of  the  a-k  curve,  i.e.. 
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if  k  is  on  the  lower  branch  of  the  a-k  curve,  then  the  solitary  impulse  uk  has 
exactly  one  unstable  mode. 

Theorem  1  (part  1).  Assume  (uko,  vkJ  is  a  solitary  impulse  solution  to  the 
system  (1)  with  k0  on  the  upper  branch  of  the  a-k  curve.  Then  for  an  arbitrary 
positive  integer  n,  the  widely  spaced  n-ple  impulse  (uk,  vk),  with  k  close  enough  to 
k0,  has  at  most  n  unstable  modes.  The  corresponding  eigenvalues  of  the  variational 
equation  approach  the  origin  as  k  -»  k0. 

Theorem  1  (part  2).  Assume  k0  is  taken  on  the  lower  branch  of  the  a-k  curve, 
and  X0  is  the  corresponding  eigenvalue  on  the  a- \  curve.  Then  the  widely  spaced 
n-ple  impulse  (uk,  vk),  with  k  close  enough  to  k0,  has  at  most  2n  -  1  and  at  least  n 
unstable  modes:  exactly  n  eigenvalues  of  the  variational  equation  approach  X0  while 
the  rest  of  them  approach  the  origin  as  k  -*  k0. 

Theorem  2  (part  1).  The  eigenvalues  close  to  the  origin,  in  both  part  1  and 
part  2  of  Theorem  1,  are  real  and  simple. 

Theorem  2  (part  2).  The  eigenvalues  close  to  X0,  in  the  part  2  of  Theorem  1 
are  real  and  simple  when  a  is  small. 

Remark  1.  By  Theorems  1  and  2  we  see  that  a  widely  spaced  n-ple  impulse 
solution  has  no  complex  unstable  mode. 


Figure  6. 
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Remark  2.  We  believe  that  the  result  of  part  2  of  Theorem  2  holds  for  all 
a  e  (0,  a„),  where  a„  is  the  knee  on  the  a-k  curve.  In  fact,  from  the  proof  of  this 
part,  we  see  that  this  is  true,  at  least,  for  a  on  the  full  interval  excluding  a  discrete 
set. 


Remark  3.  We  believe  that  the  eigenvalues  near  the  origin  correspond  to  the 
spacings,  and  those  near  \0  correspond  to  the  amplitudes.  This  idea  is  supported 
by  Feroe  [5]’s  numerical  computation.  We  recall  his  results  as  follows: 

Let  a  =  .275,  b  =  -.2,  c  =  .05,  then  the  system  (1)  has  exactly  two  solitary' 
wave  solutions:  the  larger  one  traveling  with  a  faster  velocity  than  the  smaller. 
Denote  the  faster  and  slower  velocities  by  kf  and  ks,  respectively;  see  Figure  4. 

The  system  (1)  has  countably  many  double  impulse  solutions  with  velocities 
close  to  either  kf  or  ks.  Feroe  located  four  of  them  with  velocities  kfi,  kff  kSx, 


Figure  7. 
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and  kS2  satisfying  kf  <  kfi  <  kfi,  ks  >  kS[  >  kH  and 

kfi  -  kf  ■  10“19,  kh  -  kf~  10“u, 

ks  -  kSi  *  10“15,  ks  -  kS2  ~  10’7; 

see  Figure  5. 

The  double  impulse  with  velocity  kfi  is  stable.  The  one  with  velocity  kfi  is 
unstable  with  only  one  unstable  mode,  the  eigenvalue  *  0.  This  unstable  mode  is 
of  spacing  type,  since  it  locates  at  one  of  the  pulses  (the  second  here)  and  raises 
up  the  front  of  the  pulse  while  pulling  down  the  back;  see  Figure  6.  Another 
point  of  view  concerning  this  aspect  is  that  the  unstable  mode  is  indeed  similar  to 


Figure  8. 
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the  derivative  of  the  second  pulse,  since  the  eigenvalue  is  approximately  zero. 
Therefore  this  unstable  mode  causes  a  translation  of  this  pulse. 

The  double  impulse  with  velocity  k  has  two  unstable  modes,  which  are  all  of 
amplitude  type;  see  Figure  7. 

The  double  impulse  with  velocity  kSi  has  three  unstable  modes:  one  of 
spacing  type  and  two  of  amplitude  type;  see  Figure  8. 

3.  Construction  of  Unstable  Modes 
Let 

1  +  A  —  1 

0  0 
b  c  +  X  ‘ 

k  k 

It  is  easy  to  see  that  A0  is  just  the  matrix  appearing  in  the  system  (4)  of  [15].  The 
characteristic  polynomial  is 

Q(P)  =  03  +  -  k)p2  -  (2A  +  c  +  1)0  +  J  -  j(l  +  X)(c  +  A). 

Define /J, (A),  >92(A),  /J3(A)  to  be  the  three  eigenvalues  with  >3,(0)  =  a,,  i  =  1.2.3. 
In  section  10  we  shall  prove: 

Lemma  1.  A)  >  0,  ^«/?2(A)  <  0.  ^#«j83(A)  <  0  if  and  only  if  A  e  p*. 

where  p  +  is  the  shaded  region  of  Figure  9,  with  the  left  boundary  given  by 

A  =  |  -(1  +  62  +  c)  +  ((1  +  e2  -  c)2  +  4 b)'/2  -  iOk. 
with  real  parameter  —  oo  <  0  <  so. 

Solutions  of  the  system  (6)  are  sums  of  the  exponentials  exp{  {i,z}Yr  i  = 
1, 2, 3,  where  the  vectors  Yfs  are  the  eigenvectors  of  Ax, 


H, 

(8)  y,  =  1  ,  /  =  1,2.3. 

—  P2  +  kfii  +  1  +  A 


Since  /8j  is  the  only  eigenvalue  with  non-negative  real  part,  a  bounded  solution  to 


(6)  and  (7)  on  the  interval  (-  oc,0)  must  be 


V 

0, 

(9) 

u 

V  . 

=  exp{0,z) 

1 

—  jSj2  +  Aj8|  +  1  +  X 

if  we  fix  1/(0)  =  1.  Let 


(10) 


V 

u 

v 


=  Z?,ex p{  }  K,  +  £jexp{  fi2:  }  Y2  +  Fpx p{  }  L,. 


if  zy  _,  <  z  <  zr  j  =  1,2,-  •  -  ,2n,  taking  z2„  =  +  oo.  In  view  of  the  jump  condi 
tion  (7)  we  have 


D,~  1- 


kq[ 


i 

I 

/=o 


/c/Si  +  X  +  c  C‘  exp{  -  ft, r, )  U(z,) 


K(Ol 


(ID 


kfi2  +  \  +  c  Jyl  exp {-p2z,)U(z,) 

Wi  ,_0  l“*U)l 

*A+_X  +  C  V  eXP(  —^3^  )  t/(^i) 
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where  q'  =  {?'(/).),  1  =  1-2,3.  Substituting  (8)  and  (11)  into  (10),  we  have 


U{z)  =  exp{/?,z }  -  £ 


kfis  +  \  +  c  Jy  exp {ft(z  -  z,)}  U(z,) 
k(ls  Za  K(z,)l 


if  zy_,  <  z  <  zy,  j  =  1,2,-  •  -,2m.  Therefore, 

.  r  a  i  ^  kji  +  A  +  c  ^  exp{ft(z  -  z,))u(z,) 
t/O^-ex p(^)  -  I  -  E  - - 


-exp{ft2/)+  I  *><;,, -z,)^ 


for  j  =  1,2,*  •  -  ,2n  -  1,  where 


is  (  \  ^  ^  k/i  +  \  +  c 

K\(z)-  -  L  —^7 - exp {  /?, z  }  . 


Equation  (12)  can  be  written  in  terms  of  a  matrix  as 


1 

l/U.) 


1 

(/(z,) 


exp{/?,z1} 


U(z2)  =  M„(  A)  l/(*2)  +  exp{0,z2} 


I/U2.-1)  I  exp{^,z2lI  ,} 


where  M„( A)  =  (m,  J)2nx2n  is  a  lower-triangular  matrix,  with 


m,.j  = 


=  0 


if  /  >  j. 


if  i  gy. 


Equation  (14)  implies 


1  1 
U(z  J  exp{  /?,z, } 

^(*2)  =  exp  {0,z2}  , 


1/(^2.-,) 


eXp{  PlZ2n  1  } 


where  M„(A)  =  iln  -  Mn(X)  and  I2n  is  the  In  x  2 n  identity  matrix. 
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Formula  (10)  combined  with  (11)  and  (15)  gives  the  unique  solution  to  the 
system  (6)  and  (7)  satisfying  U(  —  co)  =  0.  Therefore,  U(z)  is  an  unstable  mode  if 
and  only  if  D2n(\)  =  0. 


4.  The  Function  Tn(\,k) 

Define  Tn( A,  k)  =  Dln.  By  formulas  (11)  and  (15)  we  see  that 
T„(\,k)  =  l  +  d(\,k)X„AMn~lXi2, 
where  X„ A  and  Xn  2  are  2«-dimensional  vectors  defined  as 


*„.!  = 


1  exp  {-foi}  exp{  -  filz2 }  exp(  -&xz2n  _ , } 

l«*(0)|’  K(Zl)l  ’  l«*U2)l  '  ’  \Uk(z2n-l)\ 


Xn,2  =  (l-exp{  /3lzl } ,  exp03,z2  },•••,  exp ( 0,z2„ _ , } ) ; 
X'n  2  denotes  the  transpose  of  Xn  2,  and 

kfij  +  A  +  c 


d(\,k)  =  - 


kq{ 


The  function  Tn(-,k)  is  a  well-defined  complex  analytic  function  on  the  set 
p+.  This  function  is  the  principal  tool  in  the  study  of  the  stability  of  the  n-ple 
impulse  solutions;  A  is  an  eigenvalue  of  the  variational  equations  (6)  and  (7)  if 
and  only  if  Tn( A,  k)  =  0. 


Remark.  The  parameter  k  is  taken  so  that  the  z,,  /  =  1,2,-  ■  •  ,2n  -  1,  exist. 
Hence,  the  functions  7^( A,  k),  j  ^  n,  are  all  well  defined  when  k  is  the  traveling 
speed  of  an  n-ple  impulse  solution. 


5.  Four  Basic  Facts  Concerning  T,(\,A0) 

We  summarize  four  basic  facts  concerning  the  function  T,(A,  k0)  as  follows. 

(i)  If  k0  is  on  the  upper  branch  of  the  a-k  curve,  then  in  the  half  complex 
plane  0te  A  ^  0,  the  equation  7\(A,  k0)  =  0  has  an  unique  solution:  A  =  0. 

(ii)  If  k0  is  on  the  lower  branch  of  the  a-k  curve,  then  in  the  half  complex 
plane  91*  A  ^  0,  the  equation  7\(A,  &0)  =  0  has  two  solutions:  A  =  0  and 
A  =  A0,  where  A0  is  a  positive  number  on  the  a- A  curve;  see  Figure  3. 

(iii)  The  solution  A  =  0  of  both  (i)  and  (ii)  is  simple. 

(iv)  The  solution  A  =  A0  of  (ii)  is  simple. 

These  facts  were  given  by  Rinzel-Keller  [14],  and  Feroe  [3];  (i),  (ii)  and  (iv) 
were  obtained  numerically  and  (iii)  was  proved  algebraically  by  Rinzel-Keller 
[14],  We  assume  these  four  facts  in  the  following  discussion. 
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Lemma  3. 


6.  Behavior  of  Tn(\,k)  as  |\|  -*  oo 


lim  Tn(X,  k)  =  1 

|X|  — oo 


for  k  near  k0  and  \  e  p+. 

Proof:  The  function  Tn  is  taken  as  in  formula  (11): 

k^  +  X  +  c  2^1  exp 
"  kdi  “0  \“k(Zj)\ 

Direct  computation  with  the  aid  of  Lemma  2  (see  Section  10)  shows  that 

*0,  +  X  +  c  I  1  \  „ 


Therefore  it  is  sufficient  to  prove  that  exp{  -/J,z;  }f/(z;),  j  =*  1,2,  •  •  -  ,2w  -  1, 
are  all  bounded  as  |A|  -*  oo  for  A  e  p+.  This  fact  is  seen  from  formula  (12): 

l/(2<)  -  «P(M)  +  -  O  j^j- 

Multiplying  by  exp{  -jS^}  we  obtain 

exp{  ~~ PiZj)U(Zj)  =  1  +  £  exp{-/J,(zy z,)} 


(z  -z  xexP{~M}^(^) 


in  which 

exp{  —  Pi{zj  -  z,)}  Afx(zy  -  z-) 


+  X 

+  c  k^2  ■+■  A  ■ 

kq’2 

- 

kp  3 

4-  X 

*<?3 

- exp(ft(zy- 

=  - 

+  X 
kq{ 

+  c  ,  , 

-  +  o(exp{  - 1 

=  o| 

(vT 

■)* 

o(exp{  -  A^}) 

_  J 

\ 

-  °\ 

,  A" 

) 
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as  |X|  -»  oo  for  X  e  p+.  Thus  we  have 


1 

expf-ftzjJf/fzj) 

exp{-ftz2}t/(z2) 

i  ■ 

i 

i 

.expt-^n-l)^^-!). 

,i. 

This  makes  plain  the  fact  that  exp{-/?, Zj)U(Zj),  j  =  1,2,-  •  • , 2/i  —  1,  are  all 
bounded.  The  proof  is  finished. 

7.  Behavior  of  T„(\,k)  on  Compacts 


Lemma  4. 


lim  T„(\,k)  =  [T1(\,k0)] 


uniformly  for  X  on  arbitrary  compact  subsets  of  p+. 
Proof: 


M„_i 

*2x20-1)  M/"> 


where 


and 


M,(n)  = 


1  0 
K(z2n- 1  ~  z2n-z)  , 
\U'k(Z2»-2)\ 


*2x2(fi-l) 


Therefore, 


K(z2n-2)  K(z  2„-2  z\ ) 

K(Zo)i  l«*(zl)l 

K{zln-\)  ^{z2n-l  ~  Zl) 

K(Z0)l  K(Zl)l 


K(Z2n-2  Z2n-j) 

K(Z2„-3)l 

K(z2n  1  ~  Z2n -  3 ) 

K(Z2n-3)l 


A/.1  = 


K-i 


K<->] 


-1 
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With  this  blocking,  an  inductive  formula  of  the  function  Tn(\,  k)  is  obtained  as 
follows: 


Tn(X,k)  =  l+dXnAM-lX^2 
-  1  +  dXn.l 


eXP{  PlZ2n-l)  | 

exp{0i*2„-i}  ) 


-  1  +  dXm.XAM-}xXl.lA 

(  exp {-fil22n-2}  exp{~Plz2n-l}  \ 

\  K(Z2«-2)I  ’  l«*(*2„-l)l  / 

•[A/v(")]‘^2x2<n-i)Mn-_11^_1,2 

/  exp{-^2n_2}  \ 

\  \^k  ( Z2n  —  2 )  I  \U'k(Z2n-l)\  ) 


exp  01«2„_2}\ 
exp{^r2n_i}  / 


By  means  of  Lemma  3  of  [15],  we  see  that 


/  exp{-^z2w_2)  exp{-^z2n_i}  wi-i( exp{^22w_2)  \ 
l  K(^-2)l  ’  K(^2n-l)l  I1  1  1  l  exp {/W!}/ 


(  1  exp{ ~^1(z2n_1  -  z2n_2)} 

\  lM*(  Z2n  — 2)  I  lU*(Z2n-l)l 


r2«  -  2  )  } 


=  Ti(X,  &0)  -  1  +  o(l) 
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3S  *  k°'  UIuformly  f<>r  A  on  compact  subset  of  p\  and  that 


1,2 


I  M*j„-2)|  ’  K(*2„-i 

eMP^n-2)  o 

0  exP{fe-,}  J 


*(1  +  0(1)) 


'«-i,i 

J 


~  0  +  o(l))[r._1(A,A:)  -  llrf(^El^il£2-=i}  exP{~/W,-,]  \ 

\exP{^2,~i) 

-  r.-,a, k )r,(A, *„)  -  r._,(A, *)  _  /-.(a, *„)  +  i  +  „(I), 

Therefore,  if  we  assume  that 


k)  =  [^(A,  k0)]  n~I  +  o(l), 

then  we  see,  by  (16),  that 

T*(X’k)  =  r„_1(A,*)7’1(A,yt0)  +  0(i) 
"  fa(A.*o)]"  +  o(l). 

The  proof  is  finished. 


8.  Proof  of  Theorem  1 

case  1.  k0  is  taken  on  the  upper  branch  of  the  a-k  curve. 
StBp  1  ls  to  take  a  Positive  number  A  >  0  so  large  that 
|r„(A,k)|>A  if  |A|  £  A, 
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uniformly  for  k  near  k0.  This  can  be  done  by  means  of  Lemma  3.  In  particular, 
we  have 


Tn(\,k)*  0  if  |X|fcA. 

Step  2  is  to  take  an  arbitrarily  small  positive  number  e  >  0,  such  that  on  the 
disk 


Dt  =  {X:  |X|  s  e} 

the  function  7i(X,  k0)  has  X  =  0  as  its  unique  zero,  which  is  simple.  This  is  done 
using  fact  (i)  and  (iii)  of  Section  5. 

Step  3.  Define 


S2e  A  =  (X:  ^  0,ie  g  |X|  ^  A}. 

Then  A  is  a  compact  subset  of  p+.  By  Lemma  4,  we  have 
Tn(\,k)  =  [T1{\,k0)]"  +  o(\) 
uniformly  for  Xe!l(  A,a s  k  -*  k0.  It  follows  that 

Tn(\,k)*  0  if  XeSlliA, 
when  k  is  close  enough  to  k0. 

Step  4.  By  Cauchy’s  integral  formula  we  have 

uniformly  for  X  on  a  small  disk 

Ds=  {X:  |X|s*}, 


with  8  <  e. 

Step  5.  Applying  Rouchl’s  lemma  to  the  function  T„(\,  k)  on  the  circle 
Cs  =  {X:  |X|  =  5},  we  conclude  that  the  function  7j,(X,  k)  has  exactly  n  zeros  on 
the  disk  Ds  provided  k  is  close  enough  to  k0. 

Now  we  conclude  that  the  function  T„(\,  k)  has  at  most  n  zeros  with  positive 
real  part,  i.e.,  the  n-ple  impulse  has  at  most  n  unstable  modes  in  this  case. 

case  2.  kQ  is  taken  on  the  lower  branch  of  the  a-k  curve,  and  the  facts  (ii), 
(iii)  and  (iv)  of  Section  5  are  assumed. 
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The  proof  is  the  same  as  that  of  the  case  1,  but  with  one  more  small  disk 
centered  at  A0. 


9.  Proof  of  Theorem  2 

For  simplicity  we  prove  this  theorem  for  the  double  impulse  solutions,  and 
assume  c  —  0.  The  result  remains  true  when  c  is  small.  The  proof  for  the  n-ple 
impulse  is  by  induction. 


Proof  of  part  1 :  As  in  Section  7,  we  define 
M(l)  = 


1  0 


A/<2>  = 


1  0 


and  also 


Tau\  u\  '  ■  1  exP(~Azi) 

1  {X,k)  [k(zq)\'  K(*i)| 


[m/T1 


exP{  } 


=  ,  ■■  1  1  ^(zi)exP{  ) 

]K(z0)\  K(*i)l  K(*o)IK(*i)l 


T(2)(x  k)  =  l  +d(  1  exp{ -^1(^3  -  ^2)} 

1  (x,k)  K(*2)r  i«*(^3)i 


•iv,®ri( 


exp{^i(z3  -  z2)} 


=  1  +  d 


1  1  *(z3  -  z2)exp{  -/?i(z3  -  z2)} 


It  is  clear  that 


(17) 


lim  Tj(1)(A,  k)  =  7\(A,  k0), 

k^ka 

Urn  Tj(2)(A,  k)  =  T,{\,k 0), 

ir  —*  le~ 


uniformly  for  A  on  compact  subset  of  p+.  Therefore,  7’,(1)(7’1(2))has  a  simple  real 
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Since  we  have  assumed  that  uk  is  a  double  impulse  solution,  it  follows  from 
Section  7  of  [15]  that 


Therefore, 


S2(k)  =  0. 


lim  Arg 


I  exp{ioiz2} 

\  Pi 


(1  -  exp{  -ajZj}) 


1  +  ai^l}(eXP{^>  ”  e^P( a2^1  } ) 


=  (mod  it). 

By  (19)  and  (20),  we  see  that  (18)  holds  if  and  only  if 

So  Arg{  “S1  ~  a^Piu’(zi)) 

=  lim  Arg  {(a,  -  a2)} 

k-*k0 

*  0  (mod  it), 

but  this  is  the  case,  since 

-  a2  =  <*!  +  y  +  iu, 

with  >  0,  y  >  0,  and  w  >  0. 

Remark.  The  computation  of  (20)  is  lengthy  and  we  omit  the  details.  The 
function  S2  is  the  same  as  in  Part  1,  but  in  a  different  form. 

Step  2  is  to  refine  the  result  of  Section  7  to  get 

(21)  T2(\,  k)  =  fW(\,  k)f™(\,  k)  +  o(exp{  -j8j22}). 

The  proof  is  easy  and  therefore  we  shall  not  give  the  details. 

Step  3  is  to  fix  a  small  positive  number  5,  so  that 

(22)  k)  >  y(k)  +  e 

is  satisfied  for  all  |A|  ^  5  and  k  near  k0,  with  some  suitable  constant  t  >  0.  This 
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can  be  done  since 

*)L-o  =  ai(*)  >  y(k) 

and  both  a,  and  y  are  continuous  at  k  =  k0. 

Step  4.  It  is  seen  from  (17)  that 

lim  Xj(/c)  =  lim  X 2 ( /c )  =  0. 


In  particular,  we  have 


|M*)|<R  \\2(k)\<±8 

provided  k  is  close  enough  to  k0.  Since 

lim  z2(k)  =  +oo, 

we  may  assume  that 

|c*|exp{  —yz2 }  < 

As  k  -»  k0,  the  circles 

Q  =  {X:  |X  -  X,|  =  }|c*|exp{-yz2}} 
and 


C2  =  {X:  |X  -  X2|  =  ^|c*|exp{  -  yz2 } } 

are  disjoint,  and  on  these  circles  (22)  holds.  Applying  Rouche’s  lemma  to  formula 
(21)  on  Cj  and  C2,  respectively,  we  see  that  T2(\  k)  has  exactly  one  zero  in  each 
of  C,  and  C2.  Denote  these  zeros  by  X!  and  X2.  We  claim  that  both  X1  and  X2 
are  real.  This  is  because  T2(X,  k)  takes  real  value  if  X  is  real,  and  p+  is 
symmetric  with  respect  to  the  x-axis,  so  that  the  zeros  of  T2(  \  h)  are  either 
complex  conjugate  pairs  or  real.  From  the  proof  of  Theorem  1,  we  see  that  Xx 
and  X2  are  the  only  zeros  near  the  origin.  But  they  cannot  be  complex  conjugate 
pairs,  since  both  Xt  and  X2  are  real.  The  proof  of  part  1  of  Theorem  2  is  finished. 

Proof  of  part  2:  The  formula  (20)  of  part  1  is  replaced  by  (20)*: 

[tY2)(X0,  k)  -  7\(1)(Xo,  fc)]exp{yz2} 


(20)* 


=  2  eXP^p}  (1  -  expl-^rjjje^a,  X)|  +  o(l), 
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and  the  function  exp{y z2}S2(k)  takes  the  form 
exp{yz2},S2(A;)  =  exP{“fz2)  ^  _  exp{ -«2z, })02(a,  A:)j  +  o(l), 
where 


0j  =  otj  -  a2  -  a1exp{ajr1}  +  ajexp^Zj} 

u"(zr)  .  . 

+  )  (exp  { Qii^i }  -  expfa^j}) 

«"(*r) 


(explajz,}  -  expfa^}) 


d(\,k) 


a«i  ^^-(exptatz,}  -exp{a2z,}) 


-aa^expfajZi}  -  a2exp{a2z,})  -  -  «2) 


and 


02  =  1  + 


^expi-a^} 

u'{zi)P\ 


(exp{<Vi)  -  expfajZj}). 


If  we  can  prove  that 

(23)  lim  Arg(0!(a, &)}  #  lim  Arg(02(a,  k)}  (mod  it), 


then  the  rest  of  the  proof  is  the  same  as  that  of  part  1 .  Let  us  now  prove  (23)  for 
small  a.  Direct  computation  shows  that 


(24) 


LIM 0^0,  k)  =  2ax  1  - 
a-»0  I 


^lim(a^) 


(aj  -  a2), 


uniformly  for  k  near  k0,  where  LIM  denotes  the  set  of  limit  points.  We  shall 
prove  in  Lemma  5  of  Section  10  that  the  function  has  at  most  two  limit 
points  as  a  -*  0 : 0  and  x0,  where  x0  is  a  positive  number  depending  only  on 
parameters  b  and  c.  We  shall  also  prove  in  Lemma  5  that 
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The  computation  of  (24)  is  lengthy,  we  omit  the  details.  It  follows  from  formula 
(24)  that 

lim  Arg{0i(a,  k)}  =  Argfatj  -  a2 } 

a-»  0 

uniformly  for  k  near  k0.  Also  by  direct  computation  we  obtain 

lim  02(a,  k)  =  1  +  -  a2). 

a  —  0  Pi 

It  is  obvious  that  (23)  holds  at  this  limiting  case  and  therefore  for  the  case  of 
small  values  of  a.  The  proof  is  finished. 

Remark.  Since  both  0,  and  02  are  nice  transcendental  functions  of  a,  (23) 
should  be  true  at  least  for  a  on  the  interval  (0,  a,)  excluding  a  discrete  set.  We 
believe  that  this  discrete  set  is  empty. 

10.  Proofs  of  Lemmas  1,  2.  and  5 

Lemma  1.  Assume  |b|  <  i(l  -  c )2,  then  #«/3t(X)  >  0,  #«/l2(X)  <  0, 
9te  /)3(X)  <  0  //  X  G  p+,  where  p+  is  defined  as  in  Figure  9. 

Proof:  We  know  that 


Sttfixi 0)  =  #«aj  >  0, 

#</J2(0)  =  a2  <  0, 

■S?*03(O)  =  &4  a}  <  0. 

Therefore  in  a  neighborhood  of  zero  we  have 

8j(X)  >  0,  5f«/32(X)<0,  ^«)83(X)<0. 

This  neighborhood  can  be  extended  to  a  region.  On  the  boundary  of  this  region 
one  of  the  &(X),  i  =  1,2,3,  has  zero  real  part,  i.e.,  there  exists  a  real  number  8 
such  that 


id  -  k 


0  =  \i8I  ~  /4X|  =  det 


-1 

0 


-(1  +X) 
id 
_  b 
k 


i8  + 


1 

0 

c  +  X 
k 


—  —  [(X  +  i$k )2  +  (c  +  l  +  02)(X  +  i0X)  +  c  —  b  +  02c  , 
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Figure  10. 


which  is  equivalent  to 

X  +  iOk  =  i  -(c  +  1  +  92)  ±  ((c  +  1  +  02)2  -4 (c  -  6  +  02c))‘/2 
i.e., 

X  =  i  -(c  +  1  +  e2)  ±  ((1  +  92  -  c )  +  4b)  -  iOk. 

Taking  6  as  a  real  parameter  we  obtain  two  branches  of  continuous  curves,  as  in 
Figure  10,  which  divide  the  X-plane  into  three  parts.  In  each  part  &e  /?,(X), 
i  =  1, 2, 3,  has  a  fixed  sign,  especially  we  have 

#<0,(X)>O,  #<&(X)<0,  #*&(X)<0 

if  X  €  p+. 

Lemma  2. 

/W  =  H*  +  (*2  +  4(l+X))I/2)  +  o(l), 

&(*)  =  *(*  -  (*2  +  4(1  +  X))1/2)  +  o(l), 

MV  =  -  J  +  o{l) 

as  |X|  -*  oo  and  X  e  p+. 
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Proof: 


Let 


1)8/  -  Ax |  =  i[j32(fe/3  +  \  +  c)  -  (*/8  +  A)2 
-(l+c)(k/3  +  \)-(c-b)}. 


fix)  =  ^(jc  -  A)2(x  +  c)  -  x2  -  (1  +  c)x  -  (c  -  6), 
then  |/8/  -  /fx|  =  0  if  and  only  if  /(/r/3  +  A)  =  0.  Direct  computation  shows 

/'(*)  ilM*  ~  +  +  c)  -  2x  ~  0  +  c) 

^ x ^  ^j(jc  -  A)2(x  +  c)  -  x2  -  (1  4-  c)x  -  (c  -  b) 


Fix  |x  +  c\  -  e  for  arbitrarily  small  e  >  0,  and  let  |A|  -*  oo.  We  obtain 

1 


fix) 

fix) 


x  4  c 


4- 


°U)- 


By  Rouche’s  lemma,  we  see  that  /(x)  has  a  zero  arbitrarily  close  to  -c  as 
|A|  -»  oo.  Therefore  one  of  the  /Sj,  /82,  /?3  must  satisfy 

kfi  4-  A  =  —  c  +  o(I); 

it  is  obvious  that  fix  does  not  satisfy  this  equation  if  A  e  p+  and  (A)  -»  oo. 
Without  loss  of  generality,  we  assume  in  the  following  discussion  that 

fc/3  3  +  A=  -c  +  o(l), 

in  other  words  that 

-  i  +  o(i). 


To  estimate  /8t  and  /82,  we  take  the  characteristic  polynomial  in  another  form: 

b 


I  PI  ~  Ax  |  =  \(kfi  +  A  4-  c)  /?2  -  fc/3  -  (1  +  A)  + 


1 


=  +  A  +  c) 


(/?-r)(/?-/8-)  + 


k/8  4-  A  4-  c 
b 

kp  4-  A  4-  c 


with 


0*«i(*±(*2  +  4(l+A)),/a). 
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Again,  by  Rouche’s  lemma,  we  obtain 


£i(M  =  i(k  +  (k2  •'  4(1  +  X))1/2)  +  o(l), 


and 


&(X)  =  k(k  -  (k2  +  4(1  +  X))1/2)  +  0(1). 


Lemma  5.  (i)  The  function  ay  X  ( a )  has  at  most  two  limit  points 
and  x0,  where  x0  is  a  positive  number,  and  X(a)  is  the  a-X  curve. 

(ii)  1  ~  i)*0  *  0. 


Proof  of  Lemma  5:  Step  1.  zx/a  -» p[/ax  as  a  -*  0. 
Proof: 


1  -  ap\  =  exp{  -a^Zj } 


i.e.. 


implies 


so  that, 


log(l  -  ap{)  =  -axzx 


ap\  +  o(a)  =  axzx. 


£i 

a 


Step  2.  u'{zx)/a  -*  ax  -  p[/ax  as  a  -»  0. 
Proof: 


u'{zx)  =  oa1exp{a121}  -  £  ^exp{a,z,} 

s-l  ™ 


=  aalexp{a1z1}  -  £  ^(exp{a,Zl}  -  1), 
j-i  Fs 


as  a  -»  0  : 0 
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since  E^_i (*s/p£)  =  0.  This  gives 


=  ^expta^i} 


h  *  * 


exp  {0,0^} 


-  «! 


=  «1  - 


£i 

al  ’ 


where  6S,  s  =  1,2,3,  are  some  intermediate  values. 

Step  3.  Applying  step  1,  2  and  Lemma  2  to  the  equation 

7i(X0,  k 0)  =  0, 

we  have 

1  -  [1  +  o(fl)]exp|-2(a/X(a)  )~"[l  +  o(l)]j 

=  2(a/X(a) [l  +  o(l)]  +  *(aMa)  )2(a?  -  p[ )[l  +  o(l)] 

as  a  -»  0.  This  implies  that  aj\(a)  is  a  solution  to  following  transcendental 
equation  (with  a  small  perturbation): 

1  "  exp{”2^}  =  27TXX  +  4jcM“i  ~  S')' 

It  is  easy  to  see  that  this  equation  has  at  most  two  solutions:  x  =  0  and  x  =  x0, 
where  x0  is  a  positive  number.  Therefore  as  a  -»  0  the  function  ayX(o)  has  at 
most  two  limit  points:  0  or  x0. 

Step  4  is  to  show  that  1  -  ( p{/«i)x0  #  0. 

Proof:  If  jc0  =  ax/p'x,  then  the  transcendental  equation  of  step  3  must  be 
satisfied  by  x  =  a ix/p[,  i.e., 


so  that 


1  -  e~2  =  2  +  4 


(pO 


4aiXo  -  40,^0  +  (1  +  e  2)  =  0. 
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This  is  impossible,  since 

“h2  -  4ac”=  -16a2e~2  <  0. 
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